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ABSTRACT 

We report optical spectroscopy and high speed photometry and polarimetry of 
the INTEGRAL source IGRJ14536-5522 (=Swift J1453.4-5524). The photometry, po- 
larimetry and spectroscopy are modulated on an orbital period of 3.1564(1) hours. 
Orbital circularly polarized modulations are seen from ~ to ~ — 18 per cent, unam- 
biguously identifying IGRJ14536-5522 as a polar. The negative circular polarization 
is seen over ~ 95 per cent of the orbit, which is consistent (as viewed from Earth) 
with a single pole accretor. We estimate some of the system parameters by modeling 
the polarimetric observations. 

Some of the high speed photometric data show modulations that are consistent 
with quasi-periodic oscillations (QPOs) on the order of 5-6 minutes. Furthermore, for 
the first time, we detect the (5-6) minute QPOs in the circular polarimetry. We discuss 
the possible origins of these QPOs. In addition, we note that the source undergoes 
frequent changes between different accretion states. 

We also include details of HIPPO, a new high-speed photo-polarimeter used for 
some of our observations. This instrument is capable of high-speed, multi-filtered, si- 
multaneous all-Stokes observations. It is therefore ideal for investigating rapidly vary- 
ing astronomical sources such as magnetic Cataclysmic Variables. 

Key words: accretion, accretion discs - methods: analytical - techniques: polari- 
metric - binaries: close - novae, cataclysmic variables - X-rays: stars. 



1 INTRODUCTION 

The standard picture of a cataclysmic variable (CV) is a 
binary system consisting of a Roche lobe filling red dwarf 
(known as the secondary or the donor star) and an accret- 
ing white dwarf (the primary). CVs have orbital periods of 
typically a few hours, and mass transfer is caused by an- 
gular momentum loss - see e.g. Warner (1995) for a review 
of cataclysmic variables. Approximately 20% of the known 
CVs are magnetic cataclysmic variables (mCVs), where the 
white dwarf has a strong magnetic field (see the catalogue 



of Ritter & Kolb 2003). These are further sub-divided into 
two subtypes, namely intermediate polars (IPs) and polars, 
depending on the strength of the magnetic field of the white 
dwarf and the degree of synchronism between the white 
dwarf spin and the binary orbit- see the reviews given by 
Cropper (1990) and Patterson (1994). 

In polars, also known as AM Her systems, the white 
dwarf has a sufficiently strong magnetic field to lock the sys- 
tem into synchronous rotation and to prevent completely the 
formation of an accretion disc. Instead, the material from the 
secondary overflowing the Roche lobe initially falls towards 
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the white dwarf following a ballistic trajectory until, at some 
distance from the white dwarf, the magnetic pressure over- 
whelms the ram pressure of the ballistic stream, confining 
the flow until it eventually reaches the surface where it forms 
a hot shocked region. 

It has long been known that CVs are a significant source 
of soft (< 2 keV) and medium energy (2-10 keV) X-rays 
(e.g. Patterson et al. 1984). Recent surveys with INTE- 
GRAL show that CVs are also notable sources of hard (> 
20 keV) X-rays. A large fraction of these are made up of 
magnetic CVs and in particular IPs (see Barlow et al. 2006 
and Revnivtsev et al 2008). Interestingly, out of the five 
asynchronous polars known, two are INTEGRAL sources. 

IGRJ14536-5522 (=Swift J453.4-5524) was discovered 
as a hard X-ray source by INTEGRAL (Kuiper, Keek, 
Hermsen, Jonker & Steeghs 2006) and by Swift/BAT (Mukai 
et al. 2006). A pointed Swift/XRT observation led to the 
identification with a ROSAT all-sky survey (RASS) source 
1 RXS J145341. 1-552146, and hence to its optical identifica- 
tion (Masetti et al. 2006). Revnivtsev et al. (2008) classify 
it as an IP. 

Follow-up spectroscopy and photometry with SALT 
(Southern African Large Telescope) and with the SAAO 
(South African Astronomical Observatory) 1.9m telescope 
showed that this object belongs to a rare subtype of mag- 
netic CV: a hard X-ray bright polar or a soft intermediate 
polar (Mukai, Markwardt, Tueller, Buckley, Potter, Still et 
al. 2006). Further observations were clearly needed. 

Here we report on our optical photo-polarimetric and 
spectroscopic observations of IGRJ14536-5522. The paper 
is structured as follows: In Sect. 2 we describe the overall 
design of a new high speed photo-polarimeter used for some 
of our observations. Sect. 3 gives an account of all our ob- 
servations, followed by our analysis of the spectroscopic and 
photo-polarimetric observations in Sects. 4 and 5 respec- 
tively. We finish with a discussion and summary in Sect. 
6. 



2 THE HIGH SPEED PHOTO-POLARIMETER 
(HIPPO) 

SAAO's HIPPO was designed and built in order to re- 
place its highly successful but aging single channel equiv- 
alent, namely the UCT (University of Cape Town) photo- 
polarimeter (Cropper 1985). Its purpose is to obtain simul- 
taneous all-Stokes parameters, multi-filtered observations of 
unresolved astronomical sources. In addition, it is capable 
of high speed photo-polarimetry in order to permit inves- 
tigations of rapidly varying polarized astronomical sources. 
Of particular interest are magnetic Cataclysmic Variables 
(mCVs). This is the first refereed publication of the instru- 
ment and therefore we describe the overall instrument de- 
sign, data acquisition and reduction here. 



2.1 The optical design 

Fig. 1 shows a schematic diagram of the optical layout of the 
polarimeter. Light from the telescope first encounters a field 
lens that produces a collimated beam. Within the collimated 
beam is placed a polarizing calibration filter wheel followed 



by super-achromatic 1/4 and 1/2 wave-plates. The polar- 
izing calibration filter wheel consists of 2 linear polaroids, 
1 circular polaroid (a linear followed by a 1/4 wave plas- 
tic retarder to produce a circularly polarized beam), a Lyot 
depolarizer and an open position. These filters are used for 
calibration and efficiency measurements of the instrument 
and/or the telescope. The 1/4 and 1/2 wave-plates are also 
placed in the collimated beam in order to minimise any lat- 
eral modulation of the pupil image as the wave-plates are 
rotated. A Thompson beam-splitter then produces the ordi- 
nary and extraordinary beams. All of the above polarizing 
optics are placed before any filters or apertures, to avoid 
problems caused by metallic apertures or filters with resid- 
ual stress birefringence. 

Each beam has its own neutral density, colour filter and 
aperture wheels. The beams are focused at the aperture 
wheels by lenses at the top of each channel. Fabry lenses 
re-image the pupil onto two photo-multiplier tubes. There 
is also an eye-piece (used for initial alignment) and a dark 
slide on each channel. 

The wave-plates are contra-rotated at 10Hz and there- 
fore modulate the ordinary and extraordinary beams. The 
modulation is sufficiently rapid that errors which arise as a 
result of variable atmospheric conditions or telescope guid- 
ing modulations are much reduced. In addition, modulations 
that occur as a result of wedge shaped rotating elements, 
dirt on the rotating components or dichroism from refrac- 
tion at the element surfaces appear mostly at harmonics 
that do not affect the measurement of polarization. The re- 
maining sources of error are photon statistics, which can be 
minimised by collecting a larger number of photons, and in- 
strument/telescope polarization, which can be measured by 
observing polarized and un-polarized standard stars in com- 
bination with the calibration filters. This recipe for measur- 
ing polarimetry is based on the work of Serkowski (1974). 
Measurements of all the Stokes parameters are made simul- 
taneously from the modulated beams. As both beams are 
modulated, each provides an independent measurement of 
the polarization. Therefore, simultaneous two filter observa- 
tions are possible. Linear- and circular-only modes are also 
possible by rotating only the 1/2 or 1 /4- wave-plates respec- 
tively. 

Serkowski (1974) provides the formalism for calculating 
the Stokes parameters from O and E beams that are modu- 
lated as a result of passing through two retarders in series. In 
our case, for constantly contra-rotating 1/4 and 1/2 wave- 
plates in series, the modulated intensities are given by 

l'o = \ {j + ^Q[cos8* 8 + cos4* 4 ] + i(7[sm8* 8 - sm4* 4 ] 

-7[sin6$ 6 ]j 

' 1/1 

I E = -U+-Q[-cos8* 8 -cos4* 4 ] + 

- sm8* 8 + sm4* 4 ] +V[sm6$6]j 

for the O and E beams respectively. Where 

* 8 = 8* - 4C ( ' ) +4C ( ' ) 
<±v 6 = 6*-4C ( 5 ) +2C ,( ' ) 
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Figure 1. Optical layout of the polarimeter. Channel 1 is a copy of channel 2. 
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and 9 is the angle between the fast axis of the two wave- 
plates and referenced to the 1/4 wave-plate fast axis. C* 5 ' 
and C ( *) are the zero point constant offsets for the 1/2 
and 1/4 wave-plates respectively, and I, Q, U and V are 
the Stokes parameters. From these equations one can see 
that the linear component of the polarization is modulated 
equally at the 4th and 8th harmonics of the rotation fre- 
quency. The circular component is modulated at the 6th 
harmonic. The linear polarization is measured by adding the 
amplitudes of the 4th and 8th harmonics and, similarly by 
measuring the 6th harmonic for the circular polarization. A 
least squares algorithm is used to obtain the amplitudes and 
phases of the harmonics. Correction factors are applied to 
each harmonic in order to compensate for the fact that the 
modulated signal is made up of a finite number of bins (100) . 
Efficiency factors, to compensate for the slight wavelength 
dependence in retardance of the wave-plates and instrumen- 
tal polarization, are measured by observing polarized and 
unpolarized standard stars. 

2.2 Data acquisition and reductions 

The control and data acquisition software is written in C 
and is hosted by an industrial PC. The photometer counts 
(X2), minute and second time pulses and 1/2+1/4 wave- 
plate pulses are handled by real time C code in order to 
ensure correct and absolute timely recording of the data. 
The real time code is driven by a 1 milli-second time inter- 
val interrupt driven by a lKHz signal from the time service 
provided by the Observatory. At every interrupt, the sta- 



tus of the waveplate pulses, time pulses and the photometer 
buffers are recorded. These data are then sent to the user C 
code, where on-the-fly data reductions are performed. The 
1 milli-second data stream is also saved to disk for later off- 
line data reduction. With a data rate of lKHz, the 10Hz 
rotating wave plates are sampled 100 times per revolution. 
Therefore every 0.1 seconds a polarization measurement is 
made. Off-line data reductions permit binning of the data 
to any integer multiple of 0.1 seconds. 



3 OBSERVATIONS 

Table 1 shows a log of all the observations of IGRJ14536- 
5522. 



3.1 Optical Spectroscopy 

Spectroscopic observations of IGRJ14536-5522 were made 
during July 2007, on the 1.9-m telescope located on the 
Sutherland site of the South African Astronomical Obser- 
vatory (SAAO), using the Cassegrain spectrograph with the 
SITel CCD (1752 x 266 x 15 /xm pixels). The higher resolu- 
tion grating was chosen in order to cover the H/3 , Hell and 
H7 emission lines. The lower resolution grating was chosen 
in order to cover a broader optical wavelength range. Flat 
field spectra were obtained at the beginning and/or end of 
each night and wavelength calibration was provided by ob- 
serving a CuAr lamp approximately every 20-25 minutes. 
Spectrophotometric flux standards were also observed, al- 
lowing flux calibration. Data reductions made use of the 
standard tools available through IRAF. 
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Table 1. Table of observations. Cass Spect is the Cassegrain spectrograph, 1.9m and HIPPO are the 1.9m telescope and the High speed 
Photo-POlarimeter respectively of the South African Astronomical Observatory. OG570 and BG390 arc broad-band red and blue filters 
respectively. SALTICAM is the CCD camera of the Southern African Large Telescope. UCTCCD is the University of Cape Town CCD 
camera. RSS is the Robert Stobie Spectrograph. Observed bright and faint states are indicated by h and ' in the last column. M denotes 
observations originally published in Mukai et al. 2006. 



Date 


Telescope 


Instrument 


Spectral range / 


Resolution 


No. of 


Integration 


Datasct length 








filter 


(A) 


spectra 


times(s) 


(orbits) 


14/15 Sep 2006 


1.9m 


UCT CCD 


Clear 


- 


- 


60 


~ 0.85 hM 


15/16 Sep 2006 


1.9m 


UCT CCD 


Clear 


- 


- 


60 


~ 0.37 hM 


17/18 Sep 2006 


1.9m 


UCT CCD 


Clear 


- 


- 


60 


~ 0.72 hM 


18/19 Sep 2006 


1.9m 


UCT CCD 


Clear 


- 


- 


60 


~ 0.82 hM 


14-9 May/Jun 2006 


SALT 


RSS 


5930-7212A 


1 


31 


600 


~ 0.9 !M 


6/7 Jul 2007 


1.9m 


Cass Spect 


4000-5000A 


1 


42 


600 


~ 2.1 h 


8/9 Jul 2007 


1.9m 


Cass Spect 


4000-5000A 


1 


42 


600 


~ 2.1 h 


9/10 Jul 2007 


1.9m 


Cass Spect 


4000-5000A 


1 


39 


600 


~ 2.0 h 


10/11 Jul 2007 


1.9m 


Cass Spect 


3750-7800A 


4 


42 


500 


~ 2.0 h 


25/26 Jul 2007 


1.9m 


Cass Spect 


4000-8000A 


4 


9 


1200 


~ 1.0' 


27/28 Feb 2008 


1.9m 


HIPPO 


Unfiltered 


- 


- 


lms,0.1 


~ 1.2 h 


28/29 Feb 2008 


1.9m 


HIPPO 


OG570, BG39 


- 


- 


lms,0.1 


~ 1.2 h 


on /i t?^k onne 
zy/ 1 ret) ZUUo 


QATT 


CAT TTi^ 1 AAJ 


Ulear 






n 1 
U.l 


~ U.oO 


1/2 Mar 2008 


SALT 


SALTICAM 


Clear 






0.1 


~ 0.15 h 


1/2 Mar 2008P 


1.9m 


HIPPO 


B,I 






1ms 


~ 1.2 h 


2/3 Mar 2008 


SALT 


SALTICAM 


Clear 






0.1 


~ 0.17' 1 


2/3 Mar 2008 


1.9m 


HIPPO 


B,R 






1ms 


~ 1.2 h 


3/4 Apr 2008 


1.9m 


HIPPO 


B,I 






1ms, 0.1s 


~ 1.2 h 


5/6 Apr 2008 


1.9m 


HIPPO 


B,V,R,I 






1ms, 0.1s 


~ 1.5 h 


7/8 Apr 2008 


1.9m 


HIPPO 


Clcar,B,I 






1ms, 0.1s 


~ 1.2 h 


7/8 May 2008 


SALT 


SALTICAM 


U 






0.1 


~ 0.13' 1 


7/8 May 2008 


SALT 


SALTICAM 


B 






0.1 


~ 0.26 h 


9/10 May 2008 


SALT 


SALTICAM 


V 






0.1 


~ 0.23 h 


13/14 May 2008 


SALT 


SALTICAM 


R 






0.1 


~ 0.3' 1 


17/18 May 2008 


SALT 


SALTICAM 


B 






0.1 


~ 0.26 h 


18/19 May 2008 


SALT 


SALTICAM 


V 






0.1 


~ 0.1 h 



SALT RSS spectra were obtained when the system was 
in a lower state between May 14 and Jun 9 2006. These 
observations were presented in Mukai et al. (2006) and have 
been listed here for completeness and for comparison. 

3.2 Optical Polarimetry 

IGR J 14536-5522 was observed polarimetrically on the 
SAAO 1.9m telescope during the commissioning week of 
the HIPPO in February and March 2008 and then later in 
April 2008. The HIPPO was operated in its simultaneous 
linear and circular polarimetry and photometry mode (All- 
Stokes). White light observations (3500-9000A) were defined 
by the response of the two RCA31034A GaAs photomulti- 
plier tubes, whilst for others a broad blue band (3500-5500A) 
BG39 filter, a broad red band (5700-9000A) OG570 filter, 
or B,V,R,I filters were used. 

Several polarized (HD80558, HD111579, HD111613, 
HD147084, HD126593, HD147084, HD160529, HD298383, 
HD110984) and non-polarized (HD90156, HD100623) stan- 
dard stars (Hsu & Breger 1982 and Bastien et al. 1988) were 
observed in order to calculate the position angle offsets, in- 
strumental polarization and efficiency factors. Background 
sky polarization measurements were also taken at frequent 
intervals during the observations. Data reduction then pro- 
ceeded as outlined in the previous sections. 



3.3 Optical Photometry 

IGR J 14536-5522 was observed photometrically with the 
UCTCCD in 2006, as part of the polarimetric observations 
in 2008 and with SALTICAM in 2008. The observations arc 
not absolutely photometrically calibrated. 



4 SPECTROSCOPIC ANALYSIS 
4.1 The mean spectrum 

Fig. 2 shows the average of the spectra taken when the sys- 
tem was in a higher state (upper spectrum) and when it 
was observed to be in a lower state two weeks later dur- 
ing July of 2007 (lower spectrum). These observations were 
made with the SAAO 1.9m telescope. The higher state spec- 
trum is very typical for a polar (see e.g. QS Tel: Rosen et 
al. 1996, HU Aqr: Schwope, Mantel & Home, 1997), i.e. it 
exhibits Balmer emission lines (Ha - Hll) as well as neu- 
tral helium (He I AA 4026, 4387, 4471, 4713, 4922, 5016, 
5048, 5876, 6678, 7065A), ionised helium (He II AA 4542, 
4686, 5411). Call (3933A), blended CIII, Oil, NIB (4640A) 
and Fel (5172A). The phase-resolved low-resolution spectra 
clearly display bright and faint phases (not shown). 
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Figure 2. Average of the low resolution spectra. Upper and lower 
average spectra are made from the observations taken on 10/11 
(during a higher state) and the 25/26 (during a lower state) July 
2007, respectively. The lower state data has been magnified by 
10. 
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Figure 3. The Ha radial velocity curve of IGRJ14536-5522 folded 
on the 3.1564 hr period with arbitrary phasing. Crosses and dia- 
monds are the higher state 2007 and the lower state 2006 obser- 
vations respectively. 



ities. The higher and lower state data are distinguished by 
the crosses and diamonds respectively. Note that the higher 
state data appears to have a generally higher velocity than 
the lower state data. We attribute this to different emission 
regions contributing to the line emission when in different 
states: the irradiated face of the secondary star dominates 
the emission line profile during lower states, but the emission 
lines become multi-component during higher states (see e.g. 
UW Pic: Romero- Colmenero et al. 2003). Hence the multi- 
component nature of the emission lines have probably bi- 
ased our gamma velocity measurement due to asymmetric 
line profiles. 

In order to derive a zero point for the ephemeris we 
analysed our higher resolution spectroscopic data. Fig. 4 
shows these data, taken in July 2007, phased and folded 
on the above detected period and centered on the Hell 
4686A line. These observations show the multi-component 
nature more clearly (described in more detail below) . In par- 
ticular, a bright narrow component can be seen (indicated 
by the dashed curve), which is commonly associated with 
the irradiated face of the secondary star. We fitted multi(3)- 
Gaussian profiles to the trailed spectra and calculated the 
radial velocities of the 3 visible components as a function of 
phase. A best sinusoidal fit to the narrow component was 
used in order to derive the time of blue-red crossing and 
hence the time of inferior conjunction of the secondary star. 
Accordingly, this in turn is used to give the epoch to our 
ephemeris: 

T(HJD) = 2454290.14723(8) + 0.131517(4)S 

The number enclosed in brackets (on the epoch) is the 
formal one sigma error measurement. However we note that 
the true superior conjunction of the secondary may be offset 
by a value that is larger than the quoted error if the irradi- 
ated face of the secondary is not symmetric about the line 
of centers of the two stars. 

We measure a gamma velocity and semi-amplitude of 
—26.0(3) km/s and 80.8(4) km/s respectively for the nar- 
row component. The numbers in brackets are the one sigma 
standard deviations of a sine fit to the radial velocities. 
Henceforth, all our observations are phased on the above 
ephemeris. 



4.2 A spectroscopic ephemeris for the secondary. 

We measured the radial velocity of the Ha emission line in 
our July 2007 observations using a single Gaussian convo- 
lution method. These results were then combined with the 
same measurements made by Mukai et al (2006) of their 
June 2006 observations. A \ 2 minimization technique was 
then used in order to search for any periods. A period of 
3.1564(1) hours was detected which is consistent with that 
found by Mukai et al (2006) and recognized as the orbital 
period. The period error arises as a result of not being able 
to distinguish between aliases. 

Fig. 3 shows the Ha radial velocity curve folded on 
the 3.1564 hr period with arbitrary phasing. We measure 
a gamma velocity of —7.3(1) km/s and a semi-amplitude 
of 144.9(2) km/s. The numbers in brackets are the one 
sigma standard deviations of a sine fit to the radial veloc- 



4.3 The trailed spectra and Doppler tomograms 

Fig. 4 shows our Hell trailed spectra, the corresponding 
Doppler tomogram (Marsh & Home 1988) and the recon- 
structed trailed spectra of IGRJ14536-5522. The spectra 
were phase-folded on the ephemeris derived above, contin- 
uum subtracted and then the Doppler tomography code of 
Spruit (1998) was used. 

We have also calculated H/3 and H7 Doppler tomograms 
(not shown). These were found to be very similar to the Hell 
Doppler tomogram, but the trailed spectra do not show the 
narrower components as clearly. Therefore we present only 
the Hell observations and the corresponding interpretation. 

The Hell trailed spectra clearly display the multiple 
components that have been seen to some extent in other po- 
lars (e.g. Schwope et al. 1995 and Rosen et al. 1996). From 
Fig. 4, a narrow component (indicated by the dashed curve) 
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Figure 4. From left to right: Hell trailed spectra phase-folded on the derived ephemeris, the resulting Doppler tomogram and the 
reconstructed trailed spectra. 



is visible throughout the whole orbit and is generally recog- 
nised as emission from the irradiated face of the secondary 
star. The fact that it is visible throughout the whole or- 
bit suggests a relatively low inclination, although the radial 
velocity amplitudes are consistent with more moderate in- 
clinations (see e.g. V834 Cen: Potter et al. 2004). 

There is possibly a second narrow component that is 
most visible where it crosses and merges (indicated by diag- 
onal lines) with the narrow component from the secondary 
at phases ~ 0.25 and ~ 0.75. There also appears to be a 
broad but fainter underlining component that reaches max- 
imum blue-shift velocities of ~ —1000 km/s at phase ~ 0.5 
and it remains discernible with maximum redshift just be- 
fore phase ~ 1.0. This fits the general picture of emission 
from a ballistic accretion stream which then accelerates and 
becomes magnetically channeled before reaching the surface 
of the white dwarf (producing the broad component). 

The Hell Doppler tomogram in Fig. 4 shows the typ- 
ical features of a moderately inclined polar. In particular, 
emission is seen at the expected location of the irradiated 
face of the secondary, approximately centered on x- velocity 
~ km/s and y-velocity ~ 100 km/s. The x-velocity ap- 
pears to be slightly negative offset from 0. This can be ex- 
plained as being due to the leading edge of the secondary 
being more illuminated by the hot accretion region on the 
surface of the white dwarf, which also leads in orbital phase 
(e.g. V834 Cen: Potter et al 2004), although we cannot rule 
out that it may also be as a result of inaccuracies in the 
determination of the gamma velocity and/or the phasing of 
inferior conjunction on the secondary. The tomogram also 
shows emission at the expected location of the ballistic ac- 
cretion stream, seen to start at roughly the location of the 
secondary and curving away towards more negative x and 
y-velocities. There is also a very faint underlying component 
seen protruding from zero velocities to low negative veloci- 
ties. This has also been seen in other polars (e.g. Hu Aqr: 
Schwope, Mantel & Home 1997) and is generally explained 
as emission from accreting material where it starts to be 
threaded by the magnetic field of the white dwarf. 



5 PHOTO-POLARIMETRIC ANALYSIS 
5.1 The Photometry 

The left hand plots of Fig. 5 show the photometric observa- 
tions of IGRJ14536-5522 taken on the 5th (top 4 plots) and 
7th of April 2008 (last plot), phased on the spectroscopic 
ephemeris derived in section 4.2. The 5th of April observa- 
tions were obtained by interleaving 60 second integrations 
between the B and V filters and the R and I filters on chan- 
nels 1 and 2 of the polarimeter respectively. They are, there- 
fore, effectively simultaneous. Conditions were photometric 
during these observations. The source was in a higher state. 



5.1.1 The orbital modulation and the dip 

As can be seen from Fig. 5, the orbital modulation is single- 
humped, which can be attributed to the beaming of cy- 
clotron emission from a single accretion spot on the surface 
of the white dwarf. The orbital intensity minimum occurs 
when the emission region is most face-on. The single humped 
morphology also implies a moderate inclination for the sys- 
tem (see e.g. Ferrario & Wickramasinghe (1990) and Potter, 
Hakala & Cropper (1998)). 

The orbital intensity minimum is cut into by a narrow 
dip at approximately phase ~ 0.97. This dip is also seen 
in the polarimetry lightcurves (Fig. 5, right). Considering 
the accuracy of the phasing of the spectroscopic ephemeris 
(section 4.2), the narrow dip could be due to an eclipse of the 
white dwarf by the secondary star. However, this scenario is 
not consistent with the results of our spectroscopic analysis 
and polarimetric modeling (section 4.3 and 5.2 respectively), 
which favour moderate inclinations. Therefore, it is more 
likely that the dip arises as a result of absorption by material 
in the accretion stream and/or the accretion column directly 
above the emission region when we see it most face on (see 
e.g. Bridge et al. 2002). 
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Figure 5. The photometric and polarimctric data phased on the spectroscopic ephemeris derived in section 4.2. The upper 4 plots are 
simultaneous B,V,R and I filtered observations of the 5th of April 2008 whilst the bottom plots were made with the clear filter on the 
7th April 2008. Solid curves represent the model fit (see section 5.2). 



5.1.2 The photometric QP Ob and flickering 

A close inspection of the photometry in Fig. 5 reveals short 
period modulations throughout the orbit that are mostly 
consistent with being due to noise or flickering. However, de- 
tailed Fourier analysis of all of our photometry reveals that 
some of our data sets show significant singularly persistent 
peaks that are consistent with QPOs. 

An example data set is shown in the left (a) plot of 
Fig. 6 from the 28th February 2008. We analysed the data 
set by first splitting it into ~40 minute sections with an over- 
lap of 85 percent. Next, each section of data was normalised 
to a best fit second order polynomial before being subjected 
to Fourier analysis. The results are shown as a trailed am- 
plitude spectrum in the left (b) plot of Fig. 6. In this grey 
scale plot the darker areas correspond to larger amplitudes. 
Between phases 0.2-1.0 the amplitude spectra do not show 
any significant peaks which indicates that the variations are 
mostly flickering or noise. However, there is a significant 



dominating signal centered on 0.0032(1) Hz (5.2 minutes, 
indicated by the dashed line) between phases 1.0 and 1.3. 
This is clear evidence of a QPO. The 7 April 2008 data set 
also shows clear evidence of a QPO (Fig. 7, left (a+b) plots) 
at a similar period (0.0028(1) Hz, ~5.9(3) minutes) centered 
on the same phase range (~1.2) as well as a lower harmonic 
(0.00136(15) Hz), both indicated with dashed lines. 

In the left (c) plot of Fig. 6 we show the normalised 
photometry during the phase range that is dominated by 
the QPO. Overplotted are the least squares fit of the QPO 
frequencies. As one can see, the 28 February photometric 
QPO is very well described by the single dominant frequency 
as found in the trailed amplitude spectra. The 7 April 2008 
photometric QPO (Fig. 7, left (c) plot) shows a more vari- 
able amplitude compared to the 28 February QPO. This 
then explains why the trailed amplitude spectra shows two 
frequencies for the QPO, one being the harmonic of the 
other. The left (d) plots of Figs. 6 and 7 show the corre- 
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sponding amplitude spectra for the QPO dominated phase 
range. The QPOs are discussed further in section 6.2. 

5.2 The Polarimetry 

The right hand plots of Fig. 5 show the polarimetric (circu- 
lar) observations of IGRJ14536-5522 taken on the 5th (top 
4 plots) and 7th of April 2008 (bottom plot), phased on the 
spectroscopic ephemeris derived in section 4.2. These were 
generated from the same dataset as the photometry pre- 
sented in the left hand plots. We were unable to detect any 
linear polarization. 

The circular polarization curve is double-humped and 
negative throughout the entire orbit, which is consistent 
with emission from a single accretion region in a moderately 
inclined system. The polarization reaches maximum nega- 
tive values of ~ 15 percent, with the peaks of the humps 
occurring at orbital phases ~ 0.1 and ~ 0.7. There is also a 
narrow dip at phase ~ 0.97, where the polarization reaches 
a minimum of zero. This is coincident with the narrow dip 
seen in the photometry, thus supporting the hypothesis that 
the emission being absorbed is cyclotron radiation from the 
accretion region (section 5.1.1). The phasing of the dip im- 
plies that the location of the accretion region must be close 
to the line of centers of the two stars. In addition, there is 
a broad minimum centered on the narrow dip, which is typ- 
ically caused by the beaming of the cyclotron radiation at 
those phases where the line of sight most closely approaches 
the axis of the column (Barret & Chanmugam 1984; Wick- 
ramasinghe & Meggitt 1985). 

5.2.1 The system geometry 

We investigated the system geometry further by modeling 
the photo-polarimetry. We constructed a single arc model 
assuming accretion along dipole field lines as in Potter et 
al (1997). A magnetic field of 20MG was assumed and the 
cyclotron flux was calculated using the stratified accretion 
shock grids of Potter, Ramsay, Wu & Cropper (2002). We 
considered only the clear filter data because there does not 
appear to be any obvious wavelength dependence in the 
multi-filtered observations. Consequently, different values 
for the magnetic field strength and different shock models 
cannot be investigated. However, the results of this model 
are largely independent of our choice of the magnetic field 
strength. 

The values for the parameters describing the inclina- 
tion, dipole offset angle and location, shape and size of the 
accretion region were explored by comparing each resulting 
lightcurve to the data. The most successful set of parame- 
ters integrated the emission from an arc shaped region ex- 
tending from 170° to 230° in magnetic longitude and 10° 
from the magnetic pole with a system inclination of 50°. 
The magnetic dipole was offset by an angle of 10° from the 
white dwarf spin axis. A constant unpolarized background 
of lOOOct/s was assumed. 

We found that the model gave poor reproductions of the 
data for inclinations outside a range of 45-55°, when either 
the extent in phase or the double humped morphology of 
the circular polarization were poorly reproduced. The solid 
curves in the lower 2 plots of Fig. 5 show the results of the 



best model parameters over-plotted on the observations. The 
model has the accretion region most face on at phase ~ 0.95 
and least face on at phase ~ 0.45. 

5.2.2 The polarized QPOs and flickering 

We analysed our circular polarimetric data in the same man- 
ner as that of the photometry. The data and results for the 
28 February 2008 are shown in the right hand plots of Fig. 6. 
As can be seen from the right (b) plot, the trailed amplitude 
spectra is consistent with mostly noise and/or flickering be- 
tween phases 0.2-1.0. However a QPO centered on 0.0031(1) 
Hz (5.4(3) minutes) is clearly evident between phases 1.05- 
1.2. The QPO period is consistent within the errors to that 
of the photometric QPO. 

The right hand plots of Fig. 7 show the Fourier analy- 
sis of the circular polarimetric data taken on 7 April 2008. 
Once again, the trailed amplitude spectra displays the same 
characteristics as the photometry, i.e. a dominant frequency 
is seen to be centered on 0.0028(1) Hz (5.9(3) minutes) and 
a lower harmonic at 0.0015(1) Hz (both indicated by dashed 
lines) during the phase range 1.0-1.3. 

In the right (c) plots of Fig. 6 and Fig. 7 we show the 
normalised circular polarimetry during the phase range that 
is dominated by the QPOs. Overplotted are the least squares 
fit of the QPO frequencies. As one can see, the circularly 
polarized QPOs share the same characteristics as the pho- 
tometric QPOs, namely the 28 February circularly polarized 
QPO is very well described by the single dominant frequency 
found in the trailed amplitude spectra. Furthermore the 7 
April 2008 circularly polarized QPO shows a more variable 
amplitude compared to the 28 February QPO, thus requir- 
ing a harmonic frequency to better characterise the data. 
The QPOs are discussed further in section 6.2. 



6 DISCUSSION AND SUMMARY 

6.1 Object classification and system geometry 

Our optical spectroscopy and high speed photo-polarimetry 
of the INTEGRAL source IGRJ14536-5522 (=Swift J453.4- 
5524) unambiguously confirm its identification as a polar. 
Negative circular polarization is seen over all of the orbit, 
which is consistent with a single pole accretor at a moderate 
inclination. We estimate some of the system parameters by 
modeling the polarimetric observations. The most successful 
model integrated the emission from an arc-shaped region 
extending from 170° to 230° in magnetic longitude and 10° 
from the magnetic pole. The system inclination was found 
to be in the range of 45-55° with a magnetic dipole offset 
angle of 10°. 



6.2 The QPOs and nickering 

Our high speed photo-polarimetry shows flickering on 
minute time-scales. Furthermore, for the first time, we de- 
tect QPOs in the photometry and circular polarimetry. 

Photometric variations on the time scales of minutes are 
a common feature of many polars, from X-rays to infrared 
(e.g. Szkody & Margon 1980 and Watson, King & Williams 
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Figure 6. Left plots, a-d: The photometry, the corresponding trailed amplitude spectra, the normalised photometry for the phase range 
1.0-1.3 and its corresponding amplitude spectra respectively. The solid curve is the least squares fit using the frequency derived from the 
trailed spectra (dashed line). Right plots: as in the left plots but for the circular polarization. Data set from 28 February 2008. 



1987), and have been characterised as nickering, fluctua- 
tions, QPOs or erratic QPOs according to different authors. 
In general (although not always strictly true), QPOs de- 
scribe variations that show some coherence over a period of 
time. They mainly cluster in time-scales in the range of a 
few seconds (1-5) or minutes (4-10). The (l-5s) QPOs are of 
low amplitude and have been observed in the optical in only 
a few polars (e.g. V834 Cen, AN Uma and VV Pup). The 
larger amplitude (4-10 min) QPOs seem to be a general fea- 
ture of most polars. By eye, they can appear to show some 
sort of coherence, but fail to show any significant singularly 
persistent peaks under Fourier analysis. Instead, groupings 
of periods are seen and are often referred to as QPO-like. 
Flickering and fluctuations best describe variations that do 
not show any periodic behavior. 

There is evidence that the 4-10 minute QPOs may 
consist of the superposition of regular periodic oscillations: 
Bonnet-Bidaud, Somova & Somov (1991) report on obser- 
vations of AM Her during an intermediate brightness state 
where they find fluctuations to be nearly periodic, with a 
period increasing from 250s to 280s. They reason that they 
observed AM Her in a transition from bright-state to faint- 
state where a given accretion rate is able to excite stable 



oscillations. Alternatively, during the higher state, many 
oscillations may be present but masked by the superposi- 
tion of different modes corresponding to different accretion 
tubes. Therefore, during the transition, the accretion rate 
decreases, reducing the number of accretion tubes until fi- 
nally only one unique tube contributes to the emission. 

Ultimately, the photometric oscillations are caused by 
variations in the accretion flow. There are several theo- 
ries/models that attempt to explain the variations. King 
(1989) remarks that the (4-10 minute) periods are charac- 
teristic of a dynamical time at the photosphere of the com- 
panion. Irradiation of the region near the Li point will result 
in the formation of an ionisation front, which will tend to 
oscillate and therefore modulate the accretion. Other mod- 
els involve the capture of inhomogeneities of the accretion 
flow at the capture radius (Kuijpers and Pringle 1982) and 
accretion gate mechanisms (Patterson et al. 1981). 

The short period QPOs (l-5s) have been observed in ~ 
6 systems (e.g. V834 Cen, AN Uma and VV Pup). Langer, 
Chanmugam & Shaviv 1982) realised that the l-5s oscilla- 
tions were consistent with the cooling time-scales for white 
dwarf-radiated shock waves and thus QPOs are potentially 
powerful probes of the radiative shocks. Observations of 
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Figure 7. As in Fig 6 for the 7 April 2008 data set. 



VV Pup (Larson 1989) demonstrated conclusively that the 
QPOs arose in a region near the shock emission region on 
the white dwarf. 

From our results presented in Figs. 6 and 7 it is clear 
that IGRJ14536-5522 exhibited photometric and polarized 
QPOs (5-6 minutes) on two separate occasions. Addition- 
ally, the fact that we observe the QPOs oscillations in the 
polarimetry, unequivocally places the emission site at the 
cyclotron emitting shock region. Furthermore, the QPOs 
were seen during phases 1.0-1.3 only. We propose that the 
super-position of the oscillations arising from numerous ac- 
cretion tubes (similar to the mechanism of Bonnet-Bidaud, 
Somova & Somov 1991) applies to IGRJ14536-5522. How- 
ever, instead of a reduced accretion state leading to a re- 
duced number of accretion tubes, we suggest that for most 
of the orbit the superposition of emission from many accre- 
tion tubes leads to the observed flickering. However during 
the phase interval 1.0-1.3 we preferentially see the trailing 
edge of the accretion region as the white dwarf rotates. The 
leading part of the accretion region is effectively shielded 
from our view by its trailing edge for this phase interval. 
Consequently we observe a singular QPO corresponding to 
a singular accretion tube. 

Our observations of IGRJ14536-5522 indicate that it 



undergoes relatively frequent changes in accretion state, 
making it a good candidate for capturing it in transition 
and testing this scenario further. If QPOs are also present 
in the X-rays then IGR J 14536-5522 is an ideal source for in- 
vestigating the two dominating shock cooling mechanisms: 
i.e. bremsstrahlung and cyclotron cooling. 

Finally, we would like to note the importance of op- 
tical follow-up observations of candidate CV INTEGRAL 
sources, in particular with photometry and/or polarimetry. 
It has been noted that many CVs detected by INTEGRAL 
are IPs, both new and re-discoveries (see e.g. Barlow et 
al. 2006 and Revnivtsev et al 2008). However, as Pretorius 
(2009) points out, it is only through follow-up observations 
that unequivocal identifications can be made. A case in point 
is IGRJ14536-5522, which had been assumed to be an IP 
(see e.g. Revnivtsev et al. 2008) even though a spin period 
had not been detected. 
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